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A recently discovered inverse correlation between QSO redshift and long-term con- 
tinuum variability timescales was suggested to be the signature of microlensing on 
cosmological scales (Hawkins 1993). A general theoretical method for calculating such 
correlations is presented and applied to various lensing scenarios in the framework of 
A = Friedmann cosmologies. It is shown that the observed timescales can be strongly 
influenced by the observational limitations: the finite duration of the monitoring cam- 
paign and the finite photometric sensitivity. In most scenarios the timescales increase 
with source redshift, z s , although slower than the 1 + z s time dilation expected of in- 
trinsic variability. A decrease can be obtained for an extended source observed with 
moderate sensitivity. In this case, only lenses no further away than several hundreds 
Mpc participate in the lensing. The resulting optical depth is too small to explain the 
common long-term QSO variability unless an extremely high local lens density is as- 
sumed. These results do not support the idea that the reported inverse correlation can 
be attributed to microlensing of a uniform QSO sample by a uniform distribution of 
lenses. The possibility of using observations at various wavelengths and QSO samples 
at various positions to identify microlensing in QSO variability is also discussed. 
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The possibility that long-term QSO continuum variability 

^ carries the signature of microlensing bears upon two impor- 
tant issues: It may yield information on the amount and 
nature of compact dark matter on cosmological scales. It 
may also be used to investigate intrinsic properties of QSO, 
such as the size of the accretion disk (Rauch & Blandford 
1991) and the origin of the ubiquitous long-term variability, 
which is yet to be satisfactorily explained in terms of an 
internal QSO mechanism (Rees 1984). 

A basic problem in identifying this signature is that the 
microlensing-induced variability, if indeed it exists, may be 
masked by intrinsic variability on similar timescales. This 
problem can be addressed by a two step procedure. The 
first is to assume that the variability is due exclusively to mi- 
crolensing, choose a measurable property of the light curves, 
such as max Am, the maximal amplitude difference, and the- 
oretically predict its behavior under this assumption. The 
second step is to use the observed light curves to obtain 
upper bounds on the effect of microlensing on the property 
under study. These, in turn, can be translated into bounds 
on physical parameters such as the lens mass, Qi, the critical 
density fraction in lenses, and the size of the accretion disk. 



Such a program was carried out by Schneider (1993), who 
demonstrated by numeric simulations that microlensing in 
a flat Friedmann universe is expected to induce more high 
max Am light curves than are actually observed, and thus 
obtained upper limits on Qi in various lens mass ranges. 

The mean variability timescale is another measurable 
property of the light curve. Hawkins (1993) estimated the 
variability timescale in a sample of about 300 QSO moni- 
tored over 17 years by the width of the zero time-lag peak in 
the light curve autocorrelation function (ACF). He reported 
a systematic decrease of the typical timescales with increas- 
ing source redshift, z s , contrary to the 1 + z s time dilation 
expected for intrinsic variations in a uniform QSO sample. 
Hawkins suggested that this anti-correlation between vari- 
ability and source redshift is the signature of microlens- 
ing, but did not present detailed calculations to support 
this idea. The theoretical z s dependence of microlensing 
timescales was briefly discussed by Canizares (1982), who 
calculated the rise time for amplification by a given factor 
in the 'empty cone' cosmological model, taking into account 
corrections for amplification by more than one lens. The 
general form of the timescale's z s dependence was, however, 
not reported and the few numerical values quoted, indicat- 
ing an increase of the timescale with z s , were for a limited 
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z s range and for the go = 1/2 case only (the go dependence 
was reported to be weak). 

This work is a theoretical study of the correlation 
between the variability timescale and the source redshift 
(henceforth t-z correlations) under the hypothesis that 
long-term QSO variability over periods of years is dominated 
by microlensing. A general approximate method, indepen- 
dent of the details of the source surface brightness distribu- 
tion, is developed in section 2 for calculating the t-z correla- 
tions in A = Friedmann cosmologies. The resulting closed 
expression is validated against exact results for constant sur- 
face brightness sources. The predictions for the t-z correla- 
tions are presented and compared to the Hawkins results in 
section 3. Additional methods for detecting microlensing in 
t-z correlations are suggested. The sensitivity of the results 
to the observational procedure, the assumptions required for 
interpreting the Hawkins results as evidence for microlens- 
ing and the limitations of the method presented here, are 
discussed in section 4. 



2 CALCULATIONS 

The t-z correlations will be calculated with the assump- 
tion of a simple model whereby the typical source radius, 
£ s , and the mean transverse source velocity are indepen- 
dent of the source redshift, z s , and that the typical lens 
mass, mi, and the mean transverse lens velocity are both 
independent of the lens redshift, zt- It is assumed that the 
lenses are compact objects, that is contained within their 
Einstein radius, l e (equation A5), and that they are homo- 
geneously distributed with constant comoving density in a 
matter-dominated, A = Friedmann universe. 

An additional assumption, the shape of the amplified 
peak, is required for calculating the timescale, whether de- 
fined by an amplification threshold, by the mean FWHM or 
any other criterion. When the source is point-like, the am- 
plification is a known analytic function of the lens mass and 
impact parameter only (Paczyhski 1986). Realistic sources, 
however, are not point-like and furthermore, even if this is a 
good approximation for zt <C z s , it inevitably breaks down 
as zt approaches z s (See figure 2 below). For an extended 
source, the amplified light curve depends also on the source's 
surface brightness distribution and on the ratio between the 
source's angular size and the lens cross-section. As this ratio 
increases with zt, a larger fraction of the observed flux is un- 
lensed and the maximal amplification (ratio of lensed to un- 
lensed flux) drops to one. Rather than introduce one specific, 
but poorly motivated, surface brightness model, an approxi- 
mate geometrical method is developed below which assumes 
only the existence of a maximal cutoff on zt beyond which 
the magnification is negligible. While the lensing timescale 
thus defined cannot be directly related to an amplification 
threshold, the advantages of this method lie in its general- 
ity and in the simplicity of the geometrical interpretation 
and of the resulting expressions. Given a surface brightness 
model, an exact expression for the amplification threshold 
timescale can be derived from the approximate one by a sim- 
ple modification. It will be shown that the exact results thus 
obtained for a constant surface brightness model display the 
same general trends as those obtained by the approximate 
method. It will therefore be argued that the approximate 



method is adequate for investigating the overall properties 
of the t-z correlations and that the conclusions depend only 
weakly on the details of the source emission. 

We will assume that at any given time the observed 
QSO is undergoing amplification by no more than a single 
lens, i.e. that the optical depth for microlensing is small. In 
this case the mean timescale of an amplified peak is given 
by 

«*.)) = £ '(*<)§/*</ w 

where t(zi) is the typical lensing timescale for a lens at 
redshift zt and dP/dzt the differential probability for finding 
a lens near the line of sight, at zg. The small optical depth 
implies that the typical time between the light curve peaks 
is much larger than their width. (t(z s )) can therefore be 
interpreted as the width of the ACF zero time-lag peak, 
which measures in this case the width of a typical peak. 

We begin by considering the case where the observer 
and source are stationary and only the lens is moving. In 
the simple case of a point source, r is defined as the crossing 
time of ir£ e /2, where the ir/2 factor results from averaging 
over all impact parameters from to l e . In the case of an 
extended source (Fig. la), the Einstein radius l e is gener- 
alized to It, the distance between the two opposed points 
along the lens trajectory through the area subtended by the 
source which are no further away than l e from the area's 
boundary. The geometry describing a stationary lens and 
source and a moving observer is presented in figure lb, in 
which case r is defined as the crossing time of tt£ /2. 

In the general case the lensing timescale is determined 
by the combined motion of the observer, lens and source. 
The effective two dimensional apparent motion as measured 
by the observer is 

" (1 + zi)Dot " (1 + z.)D ot ~ t D ol ( ' 

where D a t and D os are the observer-lens and observer- 
source angular diameter distances, respectively, and v± the 
observer's transverse velocity. The contribution of v± to the 
apparent motion is obtained by noting that by the time the 
observer crosses £ , the lens appears to cross It . The ef- 
fective apparent motion, /tt e ff, can be estimated by (fieff) 1 ^ 2 
with the mean taken over the 3D velocity distribution func- 
tions of the lenses and sources and over the v± values of the 
QSO in the z s interval. Assuming that the distribution func- 
tions of the lenses and sources are independently Maxwellian 
with 2D r.m.s velocities of ut and u s , respectively, we obtain 

"•« ( "' £ * )= {((l + a)B„,) + + 

where u is related to the velocity dispersion a by u = 
•\/2(T. For a survey over a small angular area of the sky, 
(v\ Q ) = (v sina) 2 , where a is the angle between the line of 

* Eq. 2 is a generalization of eq. B9 in Kayser et al (1986) which 
applies only to zero-width caustics or to point-like lenses with 
<C i-s- (The case of setting tq = in fig. lb.) 
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Figure 1. The geometry of lensing by an extended source. The optical axis (dotted line) is the geodesic connecting the observer and 
source center at maximum amplification, a) Moving lens and source in a frame where the observer is stationary, b) Moving observer 
in a frame where the lens and source are stationary, r is the radial coordinate in a frame centered on the observer and x is the radial 
coordinate in a frame centered on the intersection of the two limiting rays (see appendix). 



sight and v a , the observer's peculiar motion relative to the 
cosmic microwave background (CMB). For an all-sky survey, 



{v\ Q ) = \v 2 . The lensing timescale is 



t(z£, z s ) 



2D iii e ff 



(4) 



The differential lensing probability is obtained by inte- 
grating over the lenses in the volume element subtended by 



dP : 



p tc\dt\2-K R\r\ 



e s + e e 



sin e d9 



(5) 



where pi is the space density of the lenses, Rt the expansion 
factor at the deflection event, rt the coordinate distance to 
the lens and 9 e , 9 S the angular diameters of the Einstein ring 
and the source, respectively. In the small angle approxima- 
tion dP is the sum of three terms 



dP = pi {^l 2 e c\dt\ + 2w^£ e £sc\dt\ + dvX 

^ U OS s 



(6) 



where dV is the differential volume subtended by the 
source^. The lensing path length It was defined so as to 
obtain the required limiting behavior of dP . In the limit of 
a point source dP reduces to the differential optical depth 



dP = pnrll c\dt\ , 



(7) 



and the definition of the lensing time coincides with that 
of an Amin = 1-34 amplification threshold. In the limit of 
a very extended source (negligible £ e ), dP reduces to the 
differential number of lenses in the volume subtended by 
the source 



dP = Pt dV . 



(8) 



Up to this point it was assumed that the source ampli- 
fication is caused by the motion of a single compact object 
across the line of sight. In the case of an extended source, 

t Eq. 6 is a generalized, differential form of the optical depth 
given in Kochanek & Lawrence (1990) 
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where 9 s /9 e ^ 1, the amplification by a single compact ob- 
ject becomes very small. Refsdal & Stabell (1991) suggest 
that even when 8 s /8 e is as high as 10, amplifications of up 
to 0.1m can be caused by a localized population of lenses, 
for example in a galactic halo. In this case the variability is 
due to the statistical Poissonian fluctuation in the number of 
lenses in the solid angle subtended by the source. The lensing 
timescale formalism presented here can address this type of 
microlensing variability without further modification, since 
the typical timescale for the Poissonian fluctuations is of the 
order of the crossing time of It and the small optical depth 
assumption remains valid since the whole lens population is 
associated with a single redshift value. 

It is convenient to transform to the dimensionless quan- 
tities: 

D 



D 



cH 



U 

Vo 



(9) 



where Ho the Hubble constant and i?H = c/Ho the Hub- 
ble radius. Using the expressions given in the appendix, all 
the relevant quantities can be expressed in terms of go, the 
deceleration parameter, zt, z s and £ s , which measures the 
deviation from the case of a point source. We obtain 



dP = ^tt e 



DoiDis 
Dos 



2b\{ 2 



Dot 

b os 



3/2 N 




1 + z t 
y/\ + 2q z e 



dz 



(10) 



where Dts the dimensionless lens-source angular diameter 
distance, and 



7T It 




(11) 



where it is assumed that the QSO survey covers a small 
angular area and that ut and u s do not depend on zt and z s . 
Finally, equations 10 and 11 are substituted into equation 1 
and the lensing timescale integral is calculated numerically. 
The relative angular size of the source 



D 



Dts Dos 



(12) 



diverges as zt approaches z s when £ s > (Fig. 2) and there- 
fore there exists a maximal lens redshift for detectable am- 
plification. This can be handled by introducing a maximal 
angular size ratio cutoff Cg, which is translated to a cut- 
off on the upper integration limit in eq. 1. Different values 
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Figure 2. the source to lens angular diameter ratio 6 s /0 e /£ s 
for go = 1/2 and z s = 1, 2, 3 and 4. The ratio diverges as zt 
approaches z s . Note that for small lens redshifts this ratio is bigger 
for the more distant QSO. 



of Cg are expected to approximate the behavior of differ- 
ent surface brightness models. Cg will be considered a free 
parameter of the calculations. 

The finite duration of the observational campaign, T, 
sets an upper limit on the detectable variation timescales. 
Variations on scales of ~ T or longer will appear as constant 
contributions to the light curve and therefore any measure of 
the light curve variability will reflect only the contribution 
of lenses from regions where r < T. This requires a modifi- 
cation of the integration limits in eq. 1, which is determined 
by introducing a maximal timescale cutoff C T . Figure 3 dis- 
plays the z dependence of r for the three cases where only 
one of the observer, lens or source is moving and for a case of 
combined motion. An upper cutoff on zt is required when the 
effective motion is due mainly to the observer and a lower 
cutoff when it is mainly due to the source motion. When 
the motion is mainly due to the lens, r has a maximum and 
this may result in the exclusion of an intermediate zt range 
from the integration interval. The lensing timescales due to 
the lens motion are significantly shorter than those due to 
the observer or source motion. In the case of combined mo- 
tion, when v± is comparable to ut and u s , the lens motion 
dominates the overall behavior of r. 

The exact amplification, A(8 s /8 e ,8t/8 e ), can be calcu- 
lated when the source surface brightness model is given ex- 
plicitly, and inverted to obtain 8t as a function of a minimal 
amplification A m i n . The amplification threshold timescale 
can be obtained from the approximate expressions above by 
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Figure 3. The lensing timescale for motion of observer only, lens 
only, source only and combined motion of all three as function of 
Z£ and z s for go = 1/2 and z s = 1 and 3. t s = 0.05, a = 125° 
and U£ = u s = 1.32. 



substituting the source size with an effective source size, now 
no longer a constant, 

l s (ze,z s ) = (^(A rain )-l)l e . (13) 



3 RESULTS 

The calculated t-z correlations depend on six unknown pa- 
rameters: go, Ho I mi, ui, u s , £ s and Cg. We will assume that 
ut = u s = 500 km/s (ut = u s = 1.32 for v = 380 km/s to- 
wards right ascension 11. 2h, declination —6° (Lubin et al 
1985)). We will assume that the QSO optical-UV contin- 
uum is emitted from the inner 20 Schwarzschild radii of the 
accretion disk, 



6 x 10 I4 M 8 cm, 



and 



0.01 



mi 



where Ms is the black hole 
_ff /(100kms _1 Mpc" 1 ) and mi 
into an upper limit on r with equation 9 



(14) 



(15) 

mass in 10 8 Mq, h = 
- mt/M@. T is translated 



0.016, 



h_ 

mi 1 yr 



(16) 



3.1 Fitting the Hawkins r— z correlations 

The QSO sample analyzed by Hawkins (1993) was taken 
from an area of 19 square degrees at right ascension 21.5h 
and declination —45° (i.e. a small angular area survey with 
a = 125°) that was monitored once a year over 17 years. 
About 1000 QSO candidates were initially selected by vari- 
ability and about 300 of those were included in the sample. 
The QSO were divided into two luminosity bins, —26 < 
M < -25 and -27 < M < -26, and into two redshift bins, 
1 < z s < 2 and 2 < z s < 3. The timescale was estimated by 
fitting the mean ACF of each of the 4 bins to the function 
6(exp( — at) — 1) + 1 where a and b are free parameters and 
t = ln(2)/a. Hawkins reports that for the low luminosity 
QSO the light curves have mean timescales of 2.3 and 1.4 
yr for the low and high redshift bins respectively. For the 
high luminosity QSO the mean timescales are 5.2 and 3.5 yr 
respectively. In both luminosity bins the timescale decreases 
by a factor of 0.7 as z s increases from 1.5 to 2.5. 

In applying the theoretical results derived above to this 
data, It will be assumed that £ s is independent of z s within 
the lm-wide bins of absolute magnitude and that the source 
is of an approximately uniform color. The implications of 
deviations from these assumptions are discussed in section 4. 

Figures 4a and 4b compare the calculated t-z corre- 
lations for the Ha wkins sample parameters with the ob- 
served ones for \Jhjm\ = 3 and 10 and a nearly point- 
like source with £ s = 0.05, corresponding to a black hole 
mass of Ms = 1.7 and 0.5, respectively. The flattening of 
the predicted lensing timescale is due to the onset of the 
maximal time cutoff C T . As \Jhjm\ increases, the typical 
lensing timescales become shorter and the cutoff comes into 
effect at higher z s . The observed t-z correlations differ from 
the calculated ones in their negative slope. The predicted 
timescales can be shortened by assuming higher lens and 
source velocities, but this results in a case similar to that 
shown in figure 4b, where the theory predicts r to be a ris- 
ing function of z s , contrary to what is observed. 

Figures 4c and 4d compare the lensing probability in the 
two cases. The onset of the maximal time cutoff excludes 
lenses at intermediate redshifts and decreases the lensing 
probability for \Jhjm\ = 3 relative to that of \Jhjm\ = 
1. It also introduces a peak in the lensing probability at a 
z s value that increases with go- The values required for 
P ~ 0.3, which is consistent with the assumption that every 
QSO is occ asiona lly undergo ing len sing, is of the order of 1 
for both \Jhjm\ = 1 and \Jhjm\ = 3. 

A decrease in r, such as seen in the observations, can 
be obtained for an extended source and a small angular size 
cuto ff. Fig ure 5a compares the Hawkins results to a case 
of v^7 mi = 10 and £ s = 1 (corresponding for example to 
mi = 0.01, h = 1 and Ms = 10). The decrease in the slope of 
the predicted timescales is caused by the small Cg = 1 cutoff. 
For go > and small enough zt values, 9 s /0 e increases with 
increasing z s , thereby causing a decrease of the maximal al- 
lowed zt, that is, the more distant QSO are lensed by nearer 
lenses (Fig. 2). In this case the observations and the pre- 
dicted timescales show a qualitative agreement for go > 1/4. 
The maximal zt values, as determined by the angular size 
cutoff, are 0.12, 0.09, and 0.06 for the cases of g = 1/4, 1/2 
and 1, respectively, implying that the microlensing objects 
are only several hundred Mpc away regardless of the source 
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Figure 4. The r—z correlations and lensing probability for the Hawkins sample (a = 125°, X = 17 yr) with £ s = 0.05 (point-like source), 
U£ = u s = 500 km/s and Cq = 10 for go — (solid line), 1/4 (dash-dotted line), 1/2 (dashed line) and 1 (dotted line). The go — 1/2 
calculated timescale (□), averaged over z s bins 1—2 and 2—3, are compared to the observed exponential fit timescale for the high (O) and 
low (A) luminosity bins, a) \J h/m\ = 3 (Mg = 1.7). b) \J h/m\ = 10 (Mg = 0.5). c, d) The specific lensing probability for (a) and (b), 
respectively. 
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Figure 5. The r—z correlations and lensing probability for the Hawkins sample (a = 125°, X = 17 yr) for an extended source with 
\JTkJra\ = 10, U£ = u s = 500 km/s. a) Approximate timescale, £ s = 1 (Mg = 10), Cq = 1. b) Amplification threshold timescale, l s = 3 
(Ms = 30), Amin = 1.1. c, d) The specific lensing probability for (a) and (b), respectively, go — (solid line), 1/4 (dash-dotted line), 
1/2 (dashed line) and 1 (dotted line). The go — 1/2 calculated timescales (□), averaged over z s bins 1—2 and 2—3, are compared to the 
observed exponential fit timescale and ACF zero crossing timescale for the high (0,+ respectively) and low (A,x) luminosity bins. 
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redshift. Consequently, the lensing probability is very small, 
P/Qi ~ 0.02 (Fig. 5c). Long term QSO variability is known 
to be very common (e.g. Hawkins & Veron 1993), and there- 
fore this value is unacceptably low even for Qi ~ 1. This 
rules out the idea of microlensed extended sources as an ex- 
planation for the Hawkins correlations, unless one assumes 
that the local density of lenses within ~ 300/t -1 Mpc of the 
solar system is at least Qi ~ 10. 

Figures 5b and 5d show the lensing timescales and prob- 
abilities as calculated for £ s = 3 and an amplification thresh- 
old of Amin = 1.1 (equation 13) with Schneider's (1993) ana- 
lytical approximation of A for a constant surface brightness 
distribution. The approximate results (figs. 5a and 5c) re- 
produce the exact ones quite well. The fact that the approxi- 
mate method requires a smaller source size is to be expected 
since it tends to over-estimate the lensing time at the lower 
zt values, before the onset of the angular size ratio cutoff, 
Cg. The overall similarity in the results, even for an extended 
source, where the approximation is crudest (for a point-like 
source the differences are negligible), demonstrates that the 
approximate method is adequate for studying the general 
trends in the t-z correlations and that these depend only 
weakly on the details of the source surface brightness. 

The slopes of the the observed t-z correlations can 
be made smaller, and closer to the calculated ones, if the 
timescales are defined to be the zero crossing point of the 
ACF. In this case, however, the timescales are increased to 
t(z s = 1.5) ~ 4, t(z s = 2.5) ~ 3 for the low luminosity bin 
and to ~ 5 and ~ 4.5 respectively for the high luminosity 
bin and are a factor of two higher than those calculated. 



3.2 Positional and wavelength dependence of 
the t — z correlations 

If the typical velocities of lenses and sources are of the same 
order of magnitude as that of the solar system relative to the 
CMB, there may be detectable differences in the timescales 
of the t-z correlations for small angular area surveys at dif- 
ferent angles to v a . Figure 6 shows that the typical timescales 
for microlensing of QSO samples in the direction of the ob- 
server's velocity are longer than for those perpendicular to 
it and that the difference becomes more marked as £ s in- 
creases. A factor of up to ~ 1.2 is obtained in the example 
of an extended source with small angular size cutoff, such 
that displays an inverse t-z correlation. The positional de- 
pendence is not large, but if it is actually discovered, it may 
provide conclusive evidence for the microlensing hypothesis. 

Another possible signature of microlensing is the dif- 
ferent response of the light curve in different bands. In the 
framework of the accretion disk model, the high and low 
energy photons are emitted from the inner and outer parts 
of the disk, respectively, and therefore have different effec- 
tive source sizes. Figure 7 shows the behavior of f for two 
bands of fixed QSO rest-frame wavelength, schematically la- 
beled as B and R, with the B band photons emitted from 
half the radius of the R band photons. When the source size 
is small, the typical timescale of the amplification peaks in 
the R band is longer than that of the B band. However, 
for an extended source with a small angular ratio cutoff, the 
trend is reversed and the R band light curve has significantly 




Figure 6. The dependence of the r— z correlations on the survey 
angle, a, for go = 1/2, U£ = u s = 1.32 and C T = oo. The case 
of a = 0° is compared to a = 90° for an extended source with 
£ s = 1 and Cg = 1 (solid line) and for a point-like source source 
with l s = 0.05 and Cg = 10 (dashed line). 

narrower amplification peaks since it is being amplified by 
nearer lenses. The mean time between the R band peaks is, 
however, expected to be longer, since the lensing probability 
is smaller. 



4 DISCUSSION AND SUMMARY 

The z s dependence of lensing-induced variability reflects an 
interplay between the lens cross-section, lens density, the 
source angular size and the combined motion of the ob- 
server, lens and source on cosmological scales. Even with the 
assumption of a uniform lens population and QSO sample, 
there remain in the problem six unknown free parameters 
(go, Ho I mi, ut, u s , £ s and Cg), some of which are partly 
degenerate (e.g. a small source and low source and lens ve- 
locities can lead to similar lensing timescales as the opposite 
combination). The following discussion will therefore be con- 
fined to the general trends in the t-z correlations and not 
to specific values of the parameters or to the bounds that 
can be set on them. 

The z s dependence of microlensing variability is dis- 
tinctly different from the 1 + z s cosmological time dilation 
of intrinsic variability. Unlike the constant gradient of the 
intrinsic, dilated timescale, the lensing timescale gradient 
decreases as z s increases, and in some cases even becomes 
negative (e.g. figures 4a and 5a). This trend appeared in all 
of the many parameter combinations that were checked. 

The observed t-z correlations can, in some situations, 
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Figure 7. The dependence of the r— z correlations on the rest- 
frame wavelength band for go = 1/2, U£ = u s = 1.32 and CV = 
oo. The B band photons are assumed to be emitted from half 
the radius of the R band photons. The timescales in the B and 
R bands are compared for an extended source of £ s = 1 and 2, 
respectively, and Cq = 1 (solid line) and for a point source of 
l s — 0.05 and 0.1, respectively, and Cq = 10 (dashed line). 

depend strongly on the observational procedure and its lim- 
itations: the finite duration of the monitoring campaign 
and the finite photometric sensitivity. When the lightcurves 
display clearly distinguishable peaks, there are two typical 
timescales: the mean peak width and the mean time between 
peaks. The results presented here apply to the peak width 
only. It is important that the method chosen to quantify 
the timescale of the lightcurve (e.g. ACF zero peak width, 
mean FWHM or power spectrum maximum) be matched to 
the definition of the timescales used in the theoretical pre- 
dictions and that these take into account the observational 
limitations. Failure to do so may lead to spurious relation- 
ships between the predictions and the data. 

The closest resemblance between the Hawkins results 
and the models investigated here is obtained for an extended 
source monitored with a moderate photometric sensitivity 
(Am ~ 0.1m) and a go > 1/4 cosmological model. The 
calculated timescales decrease with z s , as do the observed 
ones, but not as steeply (the fit was obtained by trial and 
error and it is possible that a better one exists). In this 
case the lenses that are causing the variability are within 
~ 300/t -1 Mpc of the solar system and the lensing proba- 
bility is very low unless the local lens density within this 
radius is very high (equivalent to ~ 10). The wavelength 
dependence of the t-z correlations can be used to test this 
hypothesis independently of arguments about the likelihood 
of having an atypically high local lens density. In this case 



the mean peak width in the longer wavelength is expected 
show shorter timescales, contrary to the intuitive expecta- 
tion. The sample position dependence can be also used for 
this purpose, although the decrease of the timescales with a 
may be harder to detect. 

Microlensing of point-like sources (i.e. point-like with 
respect to the typical lens mass) can be attained for 
of the order of 1 or less. The t-z correlations in this case 
are expected to rise monotonically (Fig. 4b), unlike the 
Hawkins results. This result is in agreement with that of 
Canizares (1982). The wavelength and positional depen- 
dence of the t-z correlations are not very useful for identi- 
fying this case of microlensing. The shorter wavelengths are 
predicted to vary on shorter timescales, as expected also for 
intrinsic variability, and the positional dependence may very 
well be undetectable. 

Many assumptions are involved in calculating the r- 
z correlations. The most problematic is that of a uniform 
QSO sample, as QSO are known to undergo a significant 
evolution. The QSO analyzed by Hawkins were binned by 
absolute magnitude, and this restricts the possible variance 
in the properties of objects in the same bin. Nevertheless, 
it is conceivable that QSO of similar luminosities but of dif- 
ferent cosmological epochs may have different source sizes. 
Dealing with such effects requires detailed modeling of the 
source evolution and emission mechanism. An additional z s 
dependence of £ s can be introduced by a temperature gradi- 
ent in the source, such as is expected in accretion disk mod- 
els. Observations in a given filter detect different parts of the 
source at each redshift, depending on the source radial tem- 
perature profile. This effect was studied in the specific case 
of a thin accretion disk, where the surface temperature de- 
creases as l s 3 ^ 4 (e.g. Shapiro & Teukolsky 1983) and there- 
fore i s cx (1 + 2 S ) -4 / 3 . It is found that this does not change 
the qualitative results in the case of point-like sources. For 
extended sources, the decreasing effective source size leads 
to t-z correlations that resemble those of the point sources 
shown in figure 4b. Thus, a temperature gradient increases 
the discrepancy between the theoretical prediction and the 
Hawkins results. The wavelength dependence of extended 
sources is also affected by the temperature gradient when 
the R and B bands are fixed filters in the observer's frame. 
For the same parameters as in figure 7, the R band varies 
faster than the B band for z s < 2.5, but this trend reverses 
at higher redshifts and resembles that of point-like sources. 

Another potential problem is that the definition of lens- 
ing timescale is based on the assumption of a compact lens, 
whose properties depend only on its mass and is character- 
ized by a symmetric lightcurve. This is no longer true for 
caustic networks where the lensing is a collective effect of 
many compact objects, which depends also on their distribu- 
tion and can give rise to complicated amplification behavior 
(e.g. Schneider & Weiss 1993). Finally, it is unclear whether 
comparing the calculated t-z correlations to an exponen- 
tial fit of the empirical ACF zero-lag peak is theoretically 
justified. A more cautious approach would be to fit the em- 
pirical ACF to a calculated one, which should also take into 
account overlapping amplification peaks, rather than just to 
a mean timescale. 

It is concluded that the calculated t-z correlations do 
not support the idea that the Hawkins results can be at- 
tributed to the microlensing of a uniform QSO sample by a 
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uniform distribution of lenses. The nature of the observed 
t-z correlations may reflect the details of the observational 
techniques employed no less than the underlying cosmolog- 
ical model. 



Appendix 

This appendix summarizes the identities required for ex- 
pressing the distances, timescales and probabilities of the 
microlensing in terms of the cosmological model parameters, 
the lens and source redshifts and the lens mass. 

The constant comoving lens density, the light path in- 
terval and the volume element are given by 



pi = Po(l + zf 
c\dt\ = R s 



dz 



(l + z)Vl + 2go 



(Al) 
(A2) 



and 

dV = ^ [q ° Z + j g ° ~ lK y i + 2 T - 1)]2 ™ : . (A3) 



g*(l + z)Vl + 2qoz 



dV = R 



3 z 2 (2 + z) 2 
4(1 + zf 



dz dQ, s 



qo > 



go = 0. 



(A4) 



where dfl s = w£ 2 /D 2 s is the solid angle subtended by the 
source and i?n = c/ Ho- 

The Einstein radius l e is defined as 



4Gm ( DoiDis 



c 2 D 



(A5) 



where Di 2 is the angular diameter distance between points 
1 and 2 (z\ < z 2 ), given by 

Di2 = R ^ (1 - 2g )(Gi- G 2 ) + (GiGl- G\G 2 ) ^ 



-D12 = Rh 



2g2(l + Zl )(l + Z2 ) 

(z 2 - zi)(2 + zi + z 2 ) 
2(1 + z 1 )(l + z 2 ) 2 ' 



go > 



(A7) 



go = 0. 



with G t = (1 + 2q z t ) 1/2 (Blandford & Kochanek 1987). 
The coordinate distance is given by 



-Rh [qoz + (go - l)(yi + 2g 2 - 1)] 

Ro g 2 (i + *) 

-Rh z(2 + z) 
Ro 2(1 + 2) ' 



go > 



(A8) 



(A9) 



go = 0, 



where Ro is the present-day expansion factor. 
The lens path length, It, is given by 

l t = R e r e (8 s + 6 e ) =l e + ^-ts ■ 

J-^os 



(A10) 



The observer's path length l a is given, to 1st order in 



e, by 



Rorot = (1 + zA£ e — . 



(AH) 



ro and xt are the solutions of the following equations (see 
e.g. Weinberg 1972) 

£ e = Rtxte, 
l s = R s x.e, 



ti = r \/l — kx 2 ( — xt\/\ — kr\ , 
i \/l — kx 2 — x s -\/l — krl . 



ro" 



(A12) 



where the radial coordinate systems r and x are defined in 
figure lb, R s is the scale factor at the source emission and 
k is the constant of curvature, which in a non-relativistic, 
matter dominated Friedmann universe is given by 



(A13) 



»=(£)W.,. 

In terms of the normalized radial coordinates, f = (Ro / Rr)t 
and x = (Ro/ Rn)x, the solution of equations A12 for xi is 
given by the root of 

= [4k(? e + $¥.)(¥. + + (f ~ 1) 2 ]£* - 

Wi + flfs f - 2(r1 - f 2 t )(p 2 - 1) - 
4fif s (fi + I3f s )(f s + I3fi)k]x 2 + 
[rl - f 2 ] 2 , (A14) 

with the parameter fi defined as 

(1+Z S )1 S 



P 



(1 + z t ) 



(A15) 
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